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ABSTRACT 


A discussion of various theoretical and experimental techniques for the 
calculation of long range interaction between two atomic systems at moderate 
separation is presented. Some applications of these techniques for obtaining 
gaseous properties are described. The forces between neutral molecules and 
metallic surfaces have also been discussed and numerical values of heats of 
adsorption for a number of systems have been calculated. 


INTRODUCTION 

Calculations of long range interaction between a pair of atoms or 
between a pair of molecules are needed in several fields, such as, molecular 
scattering, transport calculations in gases, atomic spectra and optics. The 
exact calculation of interaction energy requires a knowledge of the atomic 
wave functions describing the unperturbed systems in their ground states and 
the excited states. Except for simple systems like hydrogen and helium atoms, 
exact calculations cannot be made at present. Reasonably accurate estimates 
of the interaction Hamiltonian can, however, be made by making certain 
simplifying assumptions about the perturbing potential and the interacting 
systems (see ref. l). A number of approximate formulae for the. leading term 
in the interaction energy have been developed during the last 30 years. These 
formulae usually express the interaction energy in terms of such measurable 
quantities as the refractive index. Lamb shift, stopping power for fast 
charged particles and Rayleigh scattering cross section. It is the purpose 
of this report to discuss various techniques for calculating long range 
interaction between several systems. During these discussions, new techniques 
for obtaining information of the gaseous properties will be indicated. The 
subject of physically adsorbed gases on metallic surfaces will also be 
discussed and heats of adsorption calculated for several commonly used systems 


THEORY 

Attractive forces can occur between two nonpolar neutral atomic and 
molecular systems. (in this report, ^ ' we will be mainly concerned with the 
atom-atom long range attractive forces.) These forces are very weak and 
extend out to distances of the order of a few angstroms. A brief discussion 
of the theory of these interatomic forces is given below. (A more detailed 
discussion may be found in references lb and lc.) (*) 


(*)For a brief description of the atom-molecule and molecule-molecule long 
range forces, see the appendix 1. 


such 


(A) Classical Theory of Interaction 


Consider two weakly interacting spherical, atomic systems A and B with 
nuclear charges Z A and Z-q. Suppose that the systems are neutral and the 
separation R is sufficiently large 10 angstroms) to minimize charge 
overlap but not too large to bring in retardation effects resulting from the 
finite time for photon exchange. Let p A (r) be the charge element at a dis- 
tance r from the center of the atomic system A. 


V B 


The series expansion of the interaction potential, V a (R), at B due to the 
system A, in terms of its charge distribution function, p A (r), may be 
written as follows : 



V 


A 



P A (r) dx 
| R - r " 


( 1 ) 


= X v/(R), 
% 


if the charge distribution occupies a 
a limited region of space. 


where V a ^(R) 



Z 


r 


P (cos 0) dx 

J6 


for spherically symmetrical 


distribution. 


In this equation, i signifies the order of symmetry of the charge dis- 
tribution function p A (r ) and 0 is the angle between vectors r and Rv 

Thus the electrostatic potential to the second order, at B due to' A 
can be written as follows: 
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* 0 /*PaU; aT /V A 

V R »-sV< B > \A a 


P«(r) dT /*P.(r) r cos 0 dx 


R 




p A (r) r P 2 (cos 6) dx 


2R- 


Z A , M A * R o , R o ' D A • R o 
' B E 2 2E 3 


Where = /p^(r) r dv 

R = unit vector in the direction of R and 
o 


( 2 ) 


R • V R q 

is the tensor form of the quadrupole term. 

2R- 3 

The interaction Hamiltonian between systems A and B, V. (R), can now be 
written as follows: AB • 


Z A Z R 

t ab< r ) - - V + 


(Z B M A - Z A M B > - R c 
R 2 


M. 


“b - 3<M A • R o> (M B 


R o> 


E- 


V Z A D B + Z B D A> • R 


o_ + higher terms 


2R~ 


(3) 


In the first approximation, each term inequation (3) is .zero • In the second 
approximation, the first non-zero term is the dipole-dipole interaction term, 

i. e. 


M A • Mb - 3(M a ■ B o ) • B o ) 
„3 


The interaction energy in the second order approximation, which is proportional 
to the square of the interaction Hamiltonian, can now be written as follows: 
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+ 


+ 


V E > 




( 4 ) 


where C^, C^, and are constants. 

Since the dipole-dipole term decreases least rapidly with R, this will 
he the dominant term in U (R). We can thus approximate , ( * ) U AT3 (R) as 

"i n m Aij 


u ab ( e) - -J 

K 


% 

.6 


R 


AB 

6 


( 5 ) 


(A negative sign has been chosen, since the correction in the second 
approximation to the lowest eigenvalue is always negative. ) 

From equation (5), it is obvious that for two neutral systems in S 
states, the binding energy in the second approximation varies as R“6. This 
binding is called Van der Waal's force and constitutes the dominant term in 
the long range forces between neutral systems. The force constant, C AB> 
depends upon the state of the two systems and their mutual orientation. A 
classical calculation of is not feasible because of lack of precise 

knowledge of the dipole terms in equation (3). A quantum mechanical treatment 
of the interaction will now be described. 

(B) Quantum Mechanical Theory of Interaction 

Attention in the following discussion will be confined to the leading 
term in the interaction energy. The first quantum mechanical treatment of 
the dipole-dipole interactions in terms of the wavefunction describing the 
interacting systems was given by London (ref. 2a). His method was essentially 
a form of ordinary perturbation theory which will now be applied to calculate 

the interaction constant, C. D . 

AB 

According to the widely used perturbation theory (ref. 3), the perturbed 
Hamiltonian of an atomic system is written as (H q + H' ) whereas the eigenvalues 
and eigenfunctions of the perturbed system may be written as follows: 

cc oc 

E = S E ; f = 2 \ (6) 

n=0 n n=0 n 


(*)lt must, however, be pointed out that in cases involving one negative ion 
and one neutral atom (such as H - - H), R - ^ term may be of considerable 
importance in low energy scattering. 


k 


The Schroedinger equation describing the system for non-zero 
perturbation is: 


(H + H' - E) ¥ = 0 (T) 

o 

For any two atoms in nondegenerate spherically symmetrical ground states, the 

C AB 

second order change in energy, which is equivalent to — g-, can be written as 
follows : R 


E - E 
o n 


where summation S' extends over all states other than the ground state, 
For an (H - H) system; 

-M. = _ 6.499 e 2 a 5 /R 6 

For (He-He) system; 


'M . - 1.47 e 2 a 5 /R 6 

o 


The variational technique gives the following values for the long range 
interaction. 

(H - H) system (ref. 4); 

= - 6.499026 e 2 a^/R^ (ll> 

R 6 

(He-He) system (ref. 4a); 

=- 1.459 e 2 a 5 /R 6 (12! 

R 6 0 


The most widely accepted values for gg - for the two systems are as follows : 


(H - H) system: -6.499 e^a^/R^ 

o 

(He-He) system: 1.473 e 2 a^/R^ 


ref. 4b 
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Both the perturbation theory and the variational principle require a 
knowledge of the atomic wavefunctions. However, this is possible only in the case 
of simple one-electron and two-electron systems. Consequently, a number of 
alternative techniques of experimentally measuring the interaction Constant 
have been developed over the years. Some of the prominent techniques are 
described below: 


I. Methods based on oscillator strength determinations 

II. Methods based on multipole polarizability measurements 

III. Empirical methods 


Oscillator Strength Methods 


Dalgarno (ref. 5) has expressed van der Waals' interaction in terms of 
experimentally measurable oscillator strengths. His expression is as follows: 


C f^ f® 

AB _ 3 V y m n 

r 6 - 2R 6 m nfc) (e A _ jAj (]g B _ gB } (r A _ ^ B _ B } 

o m' o n o m o n ; 


where i 
follows : 


Jl n _B 

and f are electric dipole oscillator strengths ' defined as 


C - 2 /3 o£ - <) 2, r A ♦: 


o’ "I A m 1 


f = 2/3 (E - E®) If , 2 A f® 


j B n 1 


where r A and r^ are the position vectors of the interacting electrons. 

(*)ln the central field model, the oscillator strength for a single electron 
transition between (n, Z) and (n' , Z + l) states is defined as follows: 


„n', £, +. 1 _ , / _ / ^max 
n ,Z \2£ + 1 


( 1 v, » ± i K,l 


', Z + 1 


where A = Z for Z Z - 1 

= £ + 1 for & ->- & + l 


and | M 11 : 


Z + 1 


= bound-bound matrix element for the transition 


The oscillator strengths can be extended to include transitions into the 
continuum where the matrix element takes the following form: 

^ — h = bound- free matrix element 
where E = binding energy of the (n, Z) electron. 
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E, <{>, and V are the appropriate eigenvalues and eigenfunctions. One 
feature of Van der Waals forces is at once apparent from equation (13). 


Suppose 


e£ 


and E 


3 


are then positive, 


are the normal states of atoms A 

1 c 

E A > E 4 and E B > E B . Clearly — z 
mo no r6 


Hence the general theorem regarding dispersion forces : 
in their normal states attract each other. 


and B 
AB 
R 


Both f- values 


is then negative. 
Atoms (molecules) 


A theoretical calculation of oscillator strengths will require a 
knowledge of the appropriate wavefunctions and has been possible only in 
the simple hydrogen-like systems. However, oscillator strengths are related 
to a number of experimentally measurable parameters (see below) which can 
be used to calculate them. Furthermore, the existing experimental values 
in conjunction with oscillator strength Sum Rules (ref. 6) can be used to 
obtain additional oscillator strength values. 


An alternative, and more convenient, method of evaluating sums like 
equation (13) has been proposed by Bell (ref. 5a). This method does not 
require a detailed knowledge of individual f^ values, but enables one to 
compute equation (13) directly from selected oscillator strength sums. 
(Some of these sums. may be deduced directly from refractive index and 
Verdet constant measurements.) 


Some of the experimental techniques for measuring electric dipole 
oscillator strengths are described below. 

(a) Lamb Shift 

hydrogen and 
According to 
strength in 


(15) 


where k Q is the appropriate Lamb shift and other terms have been defined 
above . 

Thus a determination of k Q , should enable one to calculate f n . A direct 
measurement of Lamb shift for an atom also permits a calculation of the mean 
atomic excitation potential as indioated in (b) below. 

(b ) Stopping Power for Charged Particles 

The penetration of a beam of charged particles into matter is charac- 
terized by two mean excitation energies, I and I', associated with the 


The energy separation between the 2 Sig and 2 Pi, states in 
helium atoms is related to the relevant oscillator strengths. 
Bethe (ref. 7)j the Lamb shift is related to dipole oscillator 
the following manner. 


2L f (E - E ) Zn (E - E ) 

to k = a — a_ 

° 2 f (E - E ) 2 

n no n 
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stopping power and straggling, respectively. According to Fano (ref. 8), 
these quantities are related to the electric dipole oscillator strengths, f, 
in the following manner. 


An I = 


S' f An (E - E ) 
n n n . o 

S' f 
n n 


(16) 


An I' = 


S' f (E - E ) An (E - E ) 
n n n o n o 

S' f (E - E ) 
n n n o 


(17) 


where the summation, S', includes integration over the continum hut 
excludes n = 0 term. ?hus an experimental determination of I and I' 
should provide information on fn and vice versa* 

(c) Refractive Index of a Gas 

According to optical dispersion theory (ref, 9)9 the frequency- 
dependent refractive index, n, of an isotropic dielectric medium can be 
•written as follows : 


° 2 ■ <"o + «> 2 ■ 1 % 


“j - 


to - 1 Y. to 
J 


(18) 


where y. is the damping constant associated with frequency to. and 
J j 

N = number of electrons/cc 


to 


J 


= resonant frequency of the bound electrons 


Thus an experimental measurement of the refractive index of a gas at 
a number of frequencies gives information about its appropriate atomic 
oscillator strengths and hence the long range interaction constant, C^g. 
Bell (ref. 5a) has discussed calculation of Van der Waals force constants 
for several pairs of atoms from their oscillator strength sums, each of 
which can be derived from experimentally determined refractive index values. 


II. Multipole Polarizability Methods 

Another important measurable quantity that is often used in calculating 
Ca B is tiie polarizability of the atom. The dynamic polarizability a(v) 
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(*) 

of a physical system, having instantaneous dipole moment M, due to 
electrical waves of frequency v and electric veclor along X is given by 


a (v) = 2e 2 2 l(n|x|o)| 2 [E(n) - E(o)] 

x n [E(n) - E(o)] 2 - (hv) 2 


(19) 


where (n|x|o) is the matrix element of the transition n 0 and E(n) 
and E(o) are the energies of states n and o, respectively. 


The fact that polarizabilities are difficult to calculate need not 
disturb us for the values of a may be taken from the experiment or may be 
expressed in some other computationally convenient way (see below). London 
(ref. 26) established the following relationship between the force constant, 
C Ag, and the polarizability, a, of the interacting systems. 


"AB 


3 r A ~B 

2 h* h 


o i A a 
A B 


( 21 ) 


where I f s are the ionization potentials and a ! s are the electric 
dipole polarizabilities of the interacting atoms'. 

Approximate estimates of a(v) can he made using Pade approximants 
(ref. 12) and variational techniques (ref. 13). As an alternative to purely 
theoretical calculations, one can obtain these values from optical dispersion 
data and oscillator strength Sum Rules. Two important techniques of experi- 
mentally measuring polarizability are indicated below. 

(a) Refractive Index of Gases 

The refractive index n of a gas of number density N can be obtained 
from the Lorentz-Lorenz formula as follows: 


3 = k ttN a(v) 

n + 1 


( * ) The static polarizability, according to (19) takes the form: 

a x (0) = 2e 2 2 n j|(n|x|o)| 2 / /E(n) - E(o)j (20) 

(**) Identification of I with the ionization potential often leads to poor- 
results (ref. 10 and ref. ll). 
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For a case (gases) when n «1, we can modify (22) as follows: 


3 


(n +1) (n - 1) 
3 


= H irN a(v) 


or (n-l)=2nNa(v) 


( 23 ) 


Thus a measurement of n will enable a determination of a(v). (Also a 
measurement of n and an independent determination of a(v) is expected 
to give information on the number density and hence the pressure of the gas.) 

(b) Rayleigh Scattering Cross Section 

The scattering cross section of a gas for incident natural light is 
related to the refractive index and hence polarizability of the atom/ 
molecule of the gas in the following manner (ref. lU). 



where N = Loschmidt number. 

Thus a measurement of a(X) will permit a determination of n, which can 
be used to calculate a (v ) using equations (22) and (23). 

III. Empirical Methods 

(a) Thermal Molecular Beam Scattering .- As has been shown by Landau 
andLifshitz (ref. l4a) , the total elastic scattering cross section of low 
energy atoms/molecules from a potential of the form C g i ven by; 


o , (v) 

el. sc 


= Const 



( 25 a) 


where v is the relative velocity of the colliding atoms /molecules. Thus 

a measurement of a . at several energies can be used to calculate C nT1 

e_L • sc / $ \ AB 

for selected atom-molecule systems.' ' Several such measurements have been 

(*)Since the experimental scattering cross section is obtained from the meas- 
ured incident and scattered beam intensities, a measurement of a e ^ also 
enables a computation of the number density of the scattering gas. *f$his 
information, coupled with the temperature of the scattering gas, should enable 
calculation of the scattering gas pressure.) 
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reported for atom-atom and atom-molecule systems. The latest measured 
and the calculated values of are in excellent agreement (ref. lHb). 

The C Ag values may also be calculated from measurements of diffrential 
scattering cross section at small angles. 

/ % \ 

(b) Lennar d- Jones Potential . - An important empirical expression' 
widely used in the case of nonpolar atoms/molecules is the Lennard-Jones 
potential given below. 


U(R) = 4 e 


|(sf • (f)‘| 


(25b) 


where cr is an adjustable parameter. 
12 \ 


The 


e —rx ) term gives an approximation to the short range repulsion 


£ i) 11 


interaction whereas the ( -4 e — ^ ) is the familiar dispersion interaction. 
This expression will not be discussed any further in this report. 


Summary of Experimental Parameters Used for 
Calculating Long Range Forces 


Using the above techniques , the values of dipole polarizabilities , 
electrical dipole oscillator strength and C^g have been calculated by a 
number of workers. Some representative results are summarized in tables I, 
II, and III and IV below. 


(*)See "Molecular Theory of Gases and Liquids" by Herschfelder, Curtiss, 
and Bird (published by John Wiley and Sons, 195*0. 


11 



Table I.- Dipole Polarizabilities of Various Atoms and Ions; (*) 



(*) a Q = 0.53 x 1 (j cm. ^ 

(t) Calculated from London's formula, C„_, = j- ■ - ■ 4- ■ B a. a, 

AB l. + I A B 

A B 

In the case of hydrogen molecule, the ionization potential has been 
assumed to be same as for a hydrogen atom. If we ass ume the 
ionization potential of the H 2 molecule to be 15.6 eV, the dipole 
polarizability for the H 2 molecule is calculated to be 5.49 (a ) 3 . 
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Table II.- Oscillator Strengths for Resonance Transitions in Various 
Alkali Metals 


Alkali Metal 

f, Oscillator Strengths (Experimental) 

Li 

0.75 

Ref. 15 

Na 

O .98 


K 

0.99 

} Ref. 16 

Rb 

1.00 > 


Cs 

0.98 

Ref. 17 


Dalgarno (ref. 5) used these values of oscillator strengths to 
calculate the polarizability by using the following relation: 


f f 

a - 2 ' 2 o m S ' ^ o 

n (E - E ) 2 n (E - E ) 2 
on on 

The values thus calculated have been summarized in table I. 
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Table III (a).- Values of C AB for Hydrogen and Inert Gases 
(refs. 18 and 19) 


Nv B 

A. 

/ \(*) 
(in atomic units) 

H 

H 2 

He 

Ne 

A 

H 

6.50 

9.24 

2.82 

5.67 

20.15 

H 2 

9.24 

13.30 

4.21 

8.56 

29.60 

He 

2.82 

4.21 

1.46 

3.07 

9.88 

He 

5.67 

8.56 

3.07 

6.63 

20.63 

A 

20.15 

29.60 

9.88 

20.63 

68.08 


2 5 

(*)0ne atomic unit =f ~ 
electron volt-cnr- 


e<1 = 9*63 x 10 ^ erg-cm^ = 6.02 x 10~^ 


Table Ill(b).- Values of for Alkali Metals (refs. 5, 15a) 


Metals 

Li 

Na 

K 

Rb 

Cs 

Caa? atomic 

units 

(theoretical) ^ 

1.4 x 10 3 

1.6 x 10^ 

3.5 X 10 3 

3.8 x 10 3 

5.2 x 10 3 

C^, in atomic 
units 

(experimental ) 15 a 

— 

0.92 x 10 3 

1.59 x 10 3 

1.67 x 10 3 

3.46 x 10 3 
























Table IIl(c).- Values of C AB for* Alkali Metals and Inert Gases 
(ref. 5) 






C AB ln 

Atomic 

Units 



V 

He 

We 

A 

Kr 

Xe 

Al 

Exp 

Theory 

Exp 

Theory 

Exp 

Theory 

Theory 

Theory 

Li 

15 

22 

20 

48 

200 

200 

300 

480 

Na 

19 

25 

b2 

52 

200 

210 

330 

530 

K 

33 

33 

53 

68 

370 

280 

1+20 

680 

Rb 

27 

3b 

51 

69 

260 

280 

430 

700 

Cs 

31 

38 

59 

77 

250 

320 

480 

780 


(390) 

* 








( * ) This value has been derived from the measurements of Easterman et al. 
(ref. 20). 


Table Ill(d).- Values of C AB Calculated by Using the Combination Rule 

i. 

Cab = (CaA C Bb) 2 ( also shown are directly calculated values) 


System 

Cab (Calculated) 
in atomic units 

Cab (^ rom Table IIl(c)) 
in atomic units 

He - Li 

45.2 

22 

He - Wa 

48.3 

25 

He - K 

71.5 

33 

He - Rb 

74.5 

34 

He - Cs 

87.I 

38 


This table clearly indicates that the approximation C AB = (C/y^Gg-g) 2 
overestimates the van de Waals ' force constant for inert gas- 
alkali metal system by about a factor of 2. 
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Table IV.- Values of C AT3 for H 0 Molecule and Alkali Metals 

AJ3 d 


System 

in atomic units (ref. 19 ) 

Hg-Li 

87.6 

H 2 -Na 

96.3 

h 2 -k 

128.0 

H 2 -Rb 

132.0 

H 2 -Cs 

148.0 


Gaseous Adsorption on Metallic Surfaces 

Another important area where long range forces play an important role is 
gaseous adsorption on metallic surfaces. When neutral gaseous molecules are 
adsorbed on metallic surfaces, there are two distinctly different types of 
forces in play: The unsaturated valence forces leading to chemisorption and 

the dispersion forces leading to physical adsorption. The chemisorption sets 
in strongly at higher temperatures while physical adsorption (van der Waal’s 
adsorption) is predominant at low temperatures. We will confine our remarks 
only to physical adsorption. 

When a visitor atom is adsorbed on a metal surface, the following forces 
will come into play: (l) The dispersional force arising from the interaction 

between the instantaneous dipole of the adsorbed atom (adatom) and the metallic 
atoms; (2) the exchange force resulting from the exchange of electrons in the 
atom and the metal conduction band electrons (this is a repulsive force); 

(3) the ionic interaction^*) resulting from the electron transfer to the metal. 
Thus the total energy of interaction, E, between the adatom/molecule and the 
metal atoms can be written as follows : 


E = E (dispersion) + E' 


( 26 ) 


(*)Most of the earlier treatments of the ionic interaction ignore the 
electronic structure of the metal (refs. 21 and 22). Recently Wojciechowski 
(ref. 23), however, has included the atomic structure of the metal surface 
in the calculation of the adatom-metal ionic interaction. 
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where E' = E (exchange) + E (ionic) 

~ E (exchange) 

This is justifiable since the percentage of ionization of the adsorbing 
gases considered in the present calculations (He and H2) is not expected to 
be very high at the temperatures where physical adsorption phenomenon 
occurs. If, however, the charge transfer from an adsorbed species to the 
metallic substrate does occur (as would be the case for alkali adatoms on 
metallic surfaces, for instance), the ionic effects can be quite large. 
Bennett and Falicov (ref. 24a), being most concerned with such charge 
transfer and Pauli exclusion principle, developed a self-consistent field 
method for evaluating the energy of adsorption and the amount of charge 
transferred using a 1 electron model of the metal. They neglected all 
dispersional effects and assumed the action of an image force for the 
charge transfer problem. Good agreement with experimental results was 
obtained for K, Rb, and Cs adsorption on tungsten substrate despite the 
omission of dispersional effects. This indicates that if ionic forces do 
exist, they usually dominate intermolecular effects. This, however, is not 
expected to be the case with helium gas or H2 gas. The former is a closed 
shell atom, whereas in the latter case, a molecule of H2 may become adsorbed 
on the metal surface and then dissociate into two adsorbed H-atoms. These 
adsorbed atoms may subsequently diffuse into the metal lattice forming a 
solid solution (metal hydride). The dissociation energy (24b) of H2 
molecule is (36104 +_ 105)cm“l, thus making it a very unlikely process at 
low temperatures where physical adsorption occurs. 

Thus we have , 

E« E (dispersion) + E (exchange) ' (27) 

Pollard (ref. 24) has treated the exchange part of the mutual energy of 
an adsorbed neutral atom/molecule on a metal surface. Following a similar 
technique, we have calculated E (exchange) for hydrogen and helium gases 
adsorbed on several metallic surfaces (see appendix 2). The results are 
summarized in table V. During these calculations, it became apparent that 
the metals with larger characteristic radius, defined as the radius of a 
sphere containing one atom of the metal, have lower values of heat of 
adsorption, mainly due to increased equilibrium inter-atomic separation. 


(26a) 

(26b) 


CONCLUSIONS 

Various theoretical and experimental techniques for measuring the long 
range force constant between neutral atomic systems have been reviewed. 
Typical results for some commonly occurring gas-metal combinations have been 
summarized. During the course of discussion of various calculational tech- 
niques, new methods of measuring several gaseous properties have been 
indicated. For instance, an accurate measurement of Lamb shift in an atom 
should provide information on the stopping power of the gas atoms for charged 
particles. The measurement of refractive index of a gas can be used to 
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obtain the pressure of the gas and similarly, the total elastic scattering 
cross section in atomic beam measurements is expected to provide an accurate 
measure of the number density of the scattering gas. In the case of 
physical adsorption, the important role of exchange forces between the 
adatom electrons and the conduction band electrons has been demonstrated. 

It has been, found that the heat of physical adsorption is small for metals 
with large values of characteristic radii. 
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Table V.- Values of E(binding) = E( dispersion) + E(exchange) for Hydrogen 
and Helium Adsorbed on Several Alkali Surfaces 


No. 

System 

E( dispersion) at 
10“® cm for single 
atom (molecule )- 
atom from Table 
IIl(c) and IV 

E(disp) at 
10~8 cm 
integrated 
over entire 
metal 

(see note 2) 

E(disp) . 

min 

integrated 
over entire 
metal surface 

E( exchange) 

corresponding 

to E(disp) . 

* mm 

integrated over 
entire metal 

E(binding) over 
entire surface of 
metal 

1 

He + 

' Lithium 

(-) 13.15 ev 

(-) 0.255 ev 

(-) 0.0137 ev 

(+) 0.0044 ev 

(-) 0.0093 ev 

2 

He + 
Sodium 

(-) 15.93 ev 

(-) 0.220 ev 

(-) 0.0071 ev 

(+) 0.0022 3v 

(-) 0.0049 ev 

3 

He + 

Potassium 

(-) 18 . 04 ev 

(-) 0.1800 ev 

(-) 0.0033 ev 

(+) 0.0007 ev 

(-) 0.0026 ev 

4 

He + 

Rubidium 

(-) 20.25 ev 

(-) O.165 ev 

(-) 0.0024 ev 

(+) 0.0002 ev 

(-) 0.0022 ev 

5 

He + 
Cesium 

(-) 22.64 ev 

(-) 0.150 ev 

(-) 0.0018 ev 

(+) 0.0002 ev 

(-) 0.0016 ev 

6 

H 2 + 
Lithium 

(-) 51.47 ev 

(-) 0.870 ev 

(-) 0.0340 ev 

(+) 0.0120 ev 

1 

! 

i . 

(-) 0.0220 ev 












No. 

System 

E( dispersion) at 
10“ cm for single 
atom (molecule )- 
atom from Table 
IIl(c) and IV 

E(disp) at 
10”° cm 
integrated 
over entire 
metal 

(see note 2) 

E(disp) . 

min 

integarted 
over entire 
metal surface 

E( exchange) 

corresponding 

to E(disp) . 

min 

integrated over 
entire metal 

E( binding) over 
entire surface of 
metal 

7 

H 2 + 
Sodium 

(-) 57.80 ev 

(-) 0.770 ev 

(-) 0.0190 ev 

(+) 0.0055 ev 

(-) 0.0135 ev 

8 

H 2 + 

Potasium 

(-) 76.77 ev 

(-) 0.675 ev 

(-) 0.0100 ev 

(+) 0.0020 ev 

(-) 0.0080 ev 

9 

' H 2 + 
Rubidium 

(-) 79.20 ev 

(-) 0.630 ev 

(-) 0.0078 

(+) 0 . 00l6 ev 

(-) 0.0062 ev 

10 

H 2 + 
Cesium 

(-) 88.80 ev 

(-) 0.580 

(-) 0.0058 ev 

(+) 0.0010 ev 

(-) 0.00U8 ev 


Notes: (l) The binding energies are expressed in terms of electron volts. In order to obtain the 

Heats of Adsorption in units of calories/mole, use the following conversion factor: 

0.17 ev/molecule = 1+000 cals/mole 

(2) The values listed in the fourth column have been calculated using the formalism developed by 
Margenau and Kestner in their recent book entitled "Theory of Intermolecular Forces" 
(published by Pergamon Press, New York, 1971)* 
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APPENDIX I 


1. Atom-Molecule Long Range Forces 

The calculation of dispersion forces in this case is slightly more 
complicated than for an atom-atom case, mainly, because of the anisotrophy 
of the molecular polarizability. The dipole-dipole interaction potential 
can be written (ref. 20) in the following form: 


V E) 





+ 


A 


fli sp. (o 
2 


2 

cos 




(R) 



1 + A 


disp 2 



(A-l ) 


where P2 (cos 8) is the Legendre polynomial of the second order and 8 is 
the angle between R and the plane of molecular rotation and A(disp) is 
defined below. 


A(disp) = anisotropy factor for the dispersion potential 


11 1 
a - a 

U . 0 1 
a + 2a 


where a"^ and stand for the polarizability of the molecule parallel and 

perpendicular to the molecular axis, respectively. 

Quadrupole-dipole induction effects may further complicate the atom- 
molecule long range interaction calculations. 

2. Molecule-Molecule Long Range Forces 

In this case, in addition to the angular dependence of polarizability and 
the induction effects, one must also include electrostatic interactions 
between the two molecules. Thus the interaction energy becomes: 


Uab( r ) = U(dispersion) + U(induction) + U(electrostatic ) 


The induction effect has been calculated by Margenau (ref. 26). The 
electrostatic term may be written (ref. 27) as follows: 
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M 1 M 2 

U AT) (electrostatic ) = ~ — B(l,2) + 

AB R i 


M 1 Q 2 M 2 Q 1 
-^T-^ B(l,2) + B(2,l) 

R R 


fQ n Q ? M 0 M 0 1 

+ B(l,2) + B(l,2) + B(2,l) 

1 R ? R ? R ? J 


+ higher terms (A-2) 

where M, Q, and 0 have their usual significance and B(A,B) are the 
angular functions of the two molecules. 
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APPENDIX II 


CALCULATION OF HEATS OF PHYSICAL ADSORPTION BETWEEN 
NEUTRAL GASEOUS MOLECULES AND A METAL SURFACE 


At short distances, the interaction energy "between visitor atoms/ 
molecules and the metal atoms is comprised of two parts: (l) Attractive 

van der Waals interaction and (2) short range exchange forces arising from 
the overlap of electron clouds. In order to obtain the total potential 
energy of the system, the potential of the attractive forces must be added v / 
to the exchange term. On adding the exchange potential with the van der 
Waals potential, one obtains potential functions possessing characteristic 
minima. The depth of these minima determine the binding energy of the 
adatom and their position gives the equilibrium separation of the adatom 
from the metal surface. The values of these parameters for Hp and He gas 
in contact with several metals have been calculated by Pollard (ref. 2k). 
Using similar techniques, we have calculated equilibrium separation distances 
and the corresponding binding energies for H2 and He gas and several alkali 
metals. According to these techniques, the visiting atom is represented by 
a core of charge +e above the surface at a point r^ and an electron with 
the position coordinate r a . The solid is represented by atomic cores of 
charge Zje at points rj together with electrons at points r j . The 
wavefunction for the whole solid-atom system is taken to be: 




c + ^ ♦ c_ 1 2 


(A-3) 


where C + and C depend on the spin function of the electron in the 
visiting atom. The total energy of the system is given by: 


yV_ H l dx 

r * 

J f ¥ dx 


(A-U) 


(*)The addition of first order exchange forces to the second order dispersion 
terms may not cause great error if the exchange forces are small at the 
equilibrium distance between the adatom and the metal surface. This, indeed, 
is the situation in the cases considered here. 
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where 


H = - 


f% 2 

2m a 


Z J Z L 6 


2m ^ + 2 r 

j J J,L JL 


+■ 2 


L r AL 


2 2 

Z_ e Z T 2 

^ J V' v' e e 

J,k r Jk 


- 2 - 2 


J r Ja 


k r Ak r Aa 


, . 2 2 

+ 1 S 2 e 


2 r —»y* 

jk jk j aj 


(A-5) 


The exchange terms refer only to those electrons in the solid having the same 
spin as that of the electron in the visiting atom. The exchange interaction , 

A ( p =’ D ,J’ is g ive n "by: 


s m 


A(p .D ) = 
s m 


2 

e 

p(p s ) 

G(D m ) + ,(p s ) { 

P(P S ) - lj 

l F( V 1 


1 - F(D „> 



(A -6) 


where 


2 D 


F(D ) = i (l + Jl) 
m 2 \ a / 


m 


a 


2 D 


G( V - (t ) ; 


m 


p(p s } = 3 [ (t) Cosh (if) - Sinh (t)]/ ( if) 
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4 ip s ) = (a) [ Cosh (t )- 1 ]/^) 3 

p = radius of sphere of volume 7 - containing one metal electron of a 

o n 

given spin 

1/3 

- 2 r , where r is the characteristic radius of the metal. 

o S 

r = characteristic radius of the metal 
s 


3 gm. at. weight of metal 

^ 7r (Av. No.) (Density of metal) 


D m = e< l uilil:)rium distance of the adatom from the metal surface. Its 

numerical value is obtained from (D /p vs_ r ) values tabulated by 

m s s 

Pollard (ref. 24). 
a(He) = 0.62a 

o 

a(H ) = 0 . 95 a 
d o 

The total exchange interaction of both electrons in He and will be twice 
that for one of them when the metal contains equal numbers of both spin 
varieties and when the distance from the metal is large enough to make double 
exchange integral negligible as is the case for VDW forces. 

The various results are summarized in Tables VI to VIII. 
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Adsorbed Gas 
Atom or 
^XMolecule 

Metal Ny \ 

H 2 

He 

Li 

(-) 0.0220 ev 

(-) 0.0093 ev 

Na 

(-) 0.0135 ev 

(-) 0.001*9 ev 

K 

(-) 0.0080 ev 

(-) 0.0026 ev 

Kb 

(-) 0.0062 ev 

(-) 0.0022 ev 

Cs 

(-) 0.00^8 ev 

(-) 0.0016 ev 
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